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TABLE OF SYMBOLS AND ABBREVIATIONS 


ABBREVIATIONS: 


AIFS -ADVANCED INDIRECT FIRE SYSTEM 
ARMT -ADVANCED RAMJET MUNITIONS TECHNOLOGY 
BRL -BALLISTICS RESEARCH LABORATORY 


IODE -INTERACTIVE ORDINARY DIFFERENTIAL EQUATION 
SOLVER 


NPGS -NAVAL POSTGRADUATE SCHOOL 


SYMEGES: 
FORMULA PROGRAMS DESCRIPTION UNITS 
o< A angle of attack radians 
8 B Sideslip angle radians 
3 DE angle between radians 
body axis and 
inertial axis 
(yaw plane) 
3 GA angle between radians 
inertial axis 
and relative 
wind,yaw plane 
a LAMBDA wavelength feet 
p RO air density lbf-s2/ft4 
— = angular velocity rads/sec 
JS O angle between radians 
relative wind 
and inertial axis 
(pitch plane) 
8 THE angle between radians 
body axis and 
inertial axistyaw) 
Ca CD drag coefficient 
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Zenoeritt Gragg coefficient 
Hitt coefficient 
Slope lift coefficient 


moment coefficient 


slope moment coefficient 


damping coefficient 


centre of pressure inches 
location(from spike) 
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thrust error coefficient 
projectile diameter inches 
Drag lbf 
Drag due to wind LbE 
frequency Hz 
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transients 

mass lbm 


dynamic pressure 
relative wind 


initial velocity 
projectile 


ot 


lbf/f£t2-s2 
ft/sec 


ft/sec 


Veg 


VW (3 deg) 
VWXYCS deg) 
VWXZ (Sdeg) 
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T. INTRODUCTION 

A. ARMT OPERATIONAL CONCEPT 

Conventional artillery uses unguided projectiles to 
destroy a variety of targets including armored vehicles. 
Ballistic information and professional er peniene emer used 
to direct fire onto a target with little chance of a direct 
hit in the first few rounds. A forward observer provides 
feedback to the guns for correction until the target is 
destroyed with subsequent rounds.Generally this approach is 
limited to the 30 kilometer range and often point target 
locations are not known. 

The ARMT projectile,formerly called AIFS, changes the 
conventional concept by providing a fire and forget,boosted 
and terminally guided projectile with a highly predictable 
trajyectory.Conventional heat seekers and target identifica- 
tion systems are short range and require a blunt nose which 
has high drag. Active radar is not practical. Long ranges of 
60 kilometers are achieved by using solid fuel ramjet 
propulsion with an axisymmetric air inlet spike. An 
accurately predicted trajectory is achieved by ensuring 
thrust equals drag to produce zero axial acceleration and 
create a pseudo-vacuum trajectory during powered flight. At 
@ preset time (near ideal condition ramjet burnout) the 


inlet spike is discarded and canards are deployed for 
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terminal guidance control. A pasSive or reflected target 
energy system will be used for terminal guidance. Research 
is ongoing and sponsored by DARPA. Projectile 


characteristics are at Appendix A. 
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B. THRUST DRAG BALANCE CONCEPT 

Balanced thrust drag implies that the thrust equals the 
drag vectors during powered flight. The aim of the balance 
is to negate meteorological variations throughout the 
trajectory by reducing all force vectors to zero except 
gravity. The resulting pseudo-vacuum trajectory is highly 
predictable. The balance is accomplished by measuring axial 
acceleration and maintaining it at zero.Appendix B explores 
this concept in detail. The potential errors’ include 
transients due to wind induced reactions and measurement 
errors in the technical application of the balance in 
ft voghte. Further, any asymmetries such as aerodynamic 
misalignment’ or thrust will produce a lateral force which 
must be offset by a smali roll rate(sixth degree of 
freedom). This allows longitudinal forces to be the 
significant forces. 

The concept of thrust drag balance is best demonstrated 
in Figures 1-2 and 1-3.In Figure 1-2 the centre of gravity 
of the projectile is flying at an angle of attack (alpha) 
with a velocity vector (Vcg) relative to a fixed earth 
reference frame. The aerodynamic forces are lift(L) and 
drag(D). Assuming static static stability and neglecting 


transient dynamics, the angle of attack becomes zero, which 
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reduces lift to zero.If thrust equals drag then only the 
gravitational force exists. 

If a crosswind(Vw) is applied, theswenejgecei le cocks 
into the wind because of the aerodynamic force acting on the 
centre of pressure. Adjusting the thrust and again assuming 
static stability the aves can be offset as shown in Figure 
1-3 by applying the correct thrust.For future reference the 


component of drag normal to the plane containing Veg is Dw. 


V U 
a Vee 
u 
Lift 
Drag 
Drag Thrust 


Figure 1-2 Projectile Dynamics Thrust Equals Drag 
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Figure 1-3 Projectile Dynamics With Crosswind 


If a tailwind or headwind is applied, thrust is adjusted 
to compensate in the same manner. 

The dynamics or transient effects due to wind which 
create damped oscillatory motion have a cumulative net lift 
effect( as shown in Figure 1-4) as angle of attack isa 


proportional to lift.The area under the curve in Figure 1-4 


ie 


does not average to zero due to the damping. POrwertaeture 


reference this lift due to wind will be referred to as Lw. 





Figure 1-4 Angle Of Attack Versus Time 


Calculations integrating alpha(t) twice to demonstrate the 


above in the yaw or pitch plane are at Appendix B. 


See ©HESIS SCOPE 

It is emphasized that the scope of this thesis is to 
analyze dispersion from the pseudo-vacuum trajectory at 
terminal guidance time for the eight inch ARMT projectile.In 
colloquial terms the “basket “ from which the projectile 
Will be guided to target was determined. Further, this 


basket is based on thrust errors and wind effects for a 


2 


launch angle of 45 degrees and muzzle velocity of MACH 2.2 
at sea level. Linearized aerodynamics were used in all 


calculations. No attempt was made to establish a “footprint” - 


for projectile impact. 
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If. TRAJECTORY ANALYSIS METHODOLOGY 

meee VACUUM TRAJECTORY ANALYSIS 

A short and simple BASIC program was written to 
establish the vacuum trajectory and fuel burn time of the 
projectile in a two dimensional scenario where thrust 
exactly equaled drag. The Standard U.S. Atmosphere was used 
and the drag coefficient linearized (Reference 1]. The 
program listing and baseline data are at Appendix D. The key 
parameters and equations used are listed at Appendices A and 
D.The baseline pseudo-vacuum trajectory is shown in Figure 


1. 


Pee 1HREE-DEGREES-OF-FREEDOM MODEL 

The computer simulation was established using the 
Interactive Ordinary Die erent Equation (IODE) library 
program at NPGS. A listing of the three degrees of freedom 
equations is in the IODE specifications at Appendix E. The 
Simulation is based on wind tunnel results from Norden 
Systems (References 2 and 3] and assumes a horizontal plane 
at 30000 feet altitude.This height was chosen as significant 
fetant time is at altitude and the time of flight in this 
plane was near actual flight time. Primarily the model was 
used to study crosswind effects but resonance with the 


projectile short period was simulated. Different thrust 
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scenarios were also used. Figure Z=2 depicts the modeueem 


the X-Y plane. 


burnout time=101 secs 
burnout range=175100 feet 
burnout height=10990 feet 
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Figure 2-1 ARMT Pseudo-vacuum Trajectory 


Pit 
fii 


ny o ° 
Mn ctial) 16 + #(Cmg)@ + f(Cmg)@ = f(Cmaq) ¥ 


mx =tLsin Y — Dcos Y 
my = Lcos ¥+ Dsin Y 


J y (body) 


X (inertial ) 





ee, 


~~ x (body ) 


Figure 2-2 Three Degrees Freedom Model (Yaw Plane) 


The crosswinds were varied in intensity and duration to 
analyze crossrange dispersion. As well thrust was varied as 
a function of drag by scaling the drag. Cmq predictions were 
approximately -190 to -205 by Norden Systems (Reference 2] 
and an analysis was done to determine the effects of varying 
this damping coefficient. To further study the projectile 
characteristics, trajectories were simulated with no thrust 


fo iascs at lenuataon with varying wind profiles. Figure 2-3 , 
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shows the baseline three-degrees-of-freedom trajectory ex 


thrust equals drag with no wind. 
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Figure 2-3 Baseline Three-Degrees-Of-Freedom Trajectory 
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See LE VE-DEGREES-OF-FREEDOM MODEL 

Again the IODE program was used to solve a five degrees 
of freedom system of equations. The roll plane was neglected 
for the following reasons: 

-the program could not accept sufficient functions to 
define six degrees of freedon. 

-~the fins are expected to damp out initial spin due to 
gun rifling (900 RPM) in the first few seconds of flight. 
Sufficient spin to average projectile asymmetry errors to 
zero but not enough to significantly affect dynamics will 
remain.€ie. gyroscopic effects and Magnus effects) 

-the initial spin damping dynamics are beyond the scope 

of this thesis. 
The projectile was assumed to be axisymmetric allowing for 
ease of calculations in the yaw and pitch planes. Figures 
2-2 and 2-4 depict the yaw and pitch planes for the 
model,respectively. 

A listing of the equations is provided in the IQDE 
specifications at Appendix F. Trajectory simulations were 
run with varying wind profiles and thrust error scenarios. 
Terminal guidance time (Tc) was determined by the baseline 
five degrees of freedom trajectory shown at Figure 2-1 as 
101 seconds. This trajectory has no wind and thrust exactly 


equaled drag until discard of the air inlet spike and 
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canard deployment.All simulations assumed a 45 degree launch 
angle and Mach 2.2 muzzle velocity. Dispersion from the 
pseudo-vacuum trajectory was measured at terminal guidance 
time including the following variables: 

-X=range at 101 seconds 

-Y=crossrange at 101 seconds 

-Z=height at 101 seconds 

-Tb=burnout time 

-Vcg=speed at time 101 seconds 


Dispersion is defined as the difference between value for 


perturbed trajectories and the baseline trajectory. 
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Figure 2-4 Five-Degrees-Of-Freedom Pitch Plane Model 
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BD. SIX-DEGREES-OF-FREEDOM MODEL 

The Ballistics Research Laboratories (Aberdeen) six 
degrees of freedom model HTRAJ was obtained and adapted for 
use on the NPGS IBM 370 system. The adapted listing is 
included as Appendix G. Program compilation was achieved but 
output data was not satisfactory due primarily to input data 
syntax problems. Further,the program requires extensive 
modification to be able to vary thrust aa a function of drag 
with the ease of the previous models already described. Time 
was not available to correct these problems. Detailed 


documentation on the program HTRAJ is available through BRL. 
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TIT. TRAJECTORY RESULTS 
A. THREE DEGREES@OF PREEDCH 
1. Crosswind And Thrust Error Sensitivity 

Figures 3-1 and 3-2 show the effects of constant 
thrust errors as a percentage of drag with no wind.Recall 
that the three-degrees-of freedom trajectory is at constant 
altitude.The figures therefore represent a thrust error 
sensitivity in terms of burn time and range variations due 
to drag or thrust.The trajectories all pass nearly through 
the 218000 foot range mark at t=100 seconds, but the burnout 
ranges are Significantly different due to the varying burn 
times resulting from projectile thrust.Figure Sa 
demonstrates the obvious conclusion that the burnout time 
decreases as thrust increases. 

Figure 53-3 combines the effects of constant 
crosswinds and thrust errors to demonstrate resultant cross 
range dispersion.The constant wind velocity lines are very 
nearly linear, and the pattern is essentially a shift about 
the zero-zero point. If thrust is lesa than drag, the 
dispersion is in the direction of the wind vector.If thrust 


is greater than drag then the overcompensation produces a 


dispersion opposite to the wind vector. 
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Figure 3-1 Thrust Error Sensitivity (Range) 
Three Degrees Of Freedom 
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Figure 3-2 Thrust Error Sensitivity (Burntime) 
Three Degrees Of Freedom 
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Figure 3-3 Crosswind Sensitivity Three Degrees Of 
Freedom 
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Of interest was the dispersion due to the same wind 
profiles without thrust or mass loss. This represents an 


absolute dispersion with a given wind profile. These results 


are in Table 3.1. 


Table 3.1 No Thrust Crosswind Effects 
Three Degrees of Freedom 
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Thrust=O t=125 seconds 


vw (ft/sec) dispersion (feet) 
crossrange 
O 0 
60 1932.76 
90 293.390 
120 SI2.02 


The data clearly shows that crosswinds have a 
minimal effect in the thrusted scenarios.The dispersion in 
the thrust-equals- drag case is due to the transients (Lw) 


and is negligible as well. This is illustrated at Table 3.2. 
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Table 3.2 Dispersion Due to Net Lift 
from Transienta 


Thrust=drag 


vw (Cfeet/sec) dispersion (feet) 
crossrange 
O 0 
60 -018 
30 20127 
120 036 


Angle of attack damped to virtually zero during the 
trajectories (see Figure E--3) but the absolute maximum 
values were alao recorded.Figure 3-4 demonstrates the 
relation between angle of attack (alpha) and the wind 
velocity.Thrust errors had virtually no effect on the angle 


of attack magnitude. 
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Figure 3-4 Wind Effects On Maximum Angle Of Attack 


The magnitude of Cmq was determined to be 
approximately -200 in the flight speed ranges of the ARMT 
projectile [Reference 4]. To determine if this coefficient 
Significantly affected dispersion Cmq was varied from O to 
1.9 x 106 and the dispersion recorded. Figure 3-5 shows the 
relation between Cmq and dispersion (crossrange) with a 


constant wind of 60 feet/second and time of 125 seconds. 
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Figure 3-5 Damping Coefficient Effects 
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Careful examination of Figure 3-5 shows that dispe- 
rsion increases slightly for realistic values of Cmq and is 
therefore proportional to Cmq.Due to the small sensitivity 
of dispersion to Cmq for all practical purposes one can 
conclude dispersion is independent to Cmq.(The increased 
dispersion proportional to Cmq can be explained as the 
difference in the areas under a highly damped sinusoid 


versus a non-damped sinusoid is considerable.) 


3. Sinusoidal Wind Analysis 

To test the effects of phase angles and wave- 
lengths on the projectile dynamics a sinusoidal wind was 
input to test the model.The wind (Vw) was of the form Vo 
Sint2 PI Vo t/lambda + phase) where Vo was set at 60 
feet/second.Figures 3-6 and 3-7 show the trajectories for 
the sine and cosine winds to burnout at thrust equal to 90% 
drag.(Lambda=range and lambda=.5 range)Of note is the fact 
that the Sine wind produces a continuous dispersion while 
the cosine wind does not when thrust does not equal 
drag.This discovery led to an expanded study of different 
phase and wavelength sinusoidal wind inputs for various 
thrust scenarios.As well a simplified model was constructed 
as shown at Figure 3-8 and calculations confirmed the three 


degrees of freedom model results.Figures 3-9 and 3-10 depict 


the results of the calculations. 
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Figure 3-6 Sinusoidal Wind Effects (Sine) 
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Figure 3-7 Sinusoidal Wind Effects (Cosine) 
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Figure 3-8 Simplified Three Degrees Of Freedom Model 


The Simplified model was used to analytically 
evaluate the crossrange dispersion ‘¢y) in the cases of sine 
and cosine winds to confirm the validity of the three- 
degrees-of-freedom model.Calculations follow: 


General equations of motion in the y direction fron 


Figure 3-8 are: 


M dVy/dt =Dy which approximately equals D(Vw-Vy)/Vxo 


if Vxo=Vo=constant.Then transforming the equation we get: 


dvy/dt + D/(MVo)Vy = D/CMVo) Vw tt) 
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or in time differentiated notation: 
+ By = B Vw(t) aif D/CMVo) =B 


If a cosine wind Vw(t)=Vwo cos (2 PI Vo t/A) is applied 


and the differential equation solved we have: 


y(t)/A% =K je-Bt -cos(2 PI Vo t/xr) + 


DAsin(2 PI Vo t/”) 


2 PI m Vo2 
let K= 2 172m Vvocevwo 
DA_Vo 
1 + [4 PI DA 2 
1/2 m Voz 
Figure 3-93 graphically depicts the results as t increases.&A 
zero asymptote is observed and this conforms with model 
results. 


Now if a Sine wind Vw(t)= Vwo sin (2 Pi Vo t/A) is applied 


and the differential equation solved: 


VCty A =k DX» (l-cos(2 PI Vo t/nr)) - 
4 PI 1/2 m Vo2 


sin(2 PI VO t/~>)+4 PI(1/2 mVo2) (1-e7-Bt) 
DA 
Figure 3-10 graphically depicts these results as t 


increases.The asymptote is 


K | (4 PI 1/2 m Vo2Z) C1=e-FE) 
D> 


and this conforms to the model also. 
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Figure 3-9 Cosine Wind Analysis Results 
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Figure 3-10 Sine Wind Analysis Results 
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The constant wind case was also derived: 


Vw andthe solved differential 


VwCt) 
equation is: 


y(t)= Yet. = Vw m VosD e-Bt 


and as t increases y approaches Vw t. 


Speci#tie calculations were done to GOna i rm 
dispersion results at burnout.The calculations were based 


on: 


y=dy/dx X = 1/Vo dy/dt X 

Vwo=60 ft/sec 

Vo=2190 ft/sec 

‘A =304000 ft 

t= A/Vo 

Daiy2 9 Vot A Cd (=226 lbf) 
Thrust=.9D therefore D=22.6 I1bf 


M=225 lbm 


Taking the time derivative of the cosine wind y(t) 


dispersion formula derived earlier we get 


dy/dt=K|-Be-~Bt + 2 PI Vo sin(2 PI Vo t/a) + 


A 
DA oe oes 2 i/o EA 
4 PI 1/2 mVo2 a 


at t= A/Vo the equation becomes: 
dy/dt = K(-Be7~BA/VO + B) 


and substituting the values into the approximationn 
FORTY;  y=O 2c. 
Now, performing the same operations on the sine 


wind dispersion equation : 


dy/dt=K D‘A . 2 PI Vo sint2 PEvanr > 
4 PI 1/2 mVo2 ~ 


-2 PI Vo cos(2_PI Vo t) +(4 PI 1/2 m Vo), 
us n DA 


(Be~Bt) 
and at t=MWVo 
dy/dt =K{ -2_ PI Vo -4 PI 1/2 m Vo2#-eBe~Bt 
~ DX” 
and y=20 £ft approximately. 


Finally, taking the constant wind dispersion 


formula and taking its derivative 
dy/dt = Vw- Vw e7Bt 


and y= 10 £t approximately. 
These results are consistent with the three-degrees-of- 
freedom model considering the initial approximations. 
To study the short period resonance phenomenom the 


phase and wavelength wind variations were documented in 


detail for a condition of thrust equal to 90% drag.Figures 
3-11 to 3-13 show the crossrange dispersion effects fron 
wavelengths of S000 ¢€ lambda = .016 range) to 
303000(lambda=range) feet andthe phase angle variations 
from O to 2 PI.At approximately 9850 feet the dispersion 
characteristics were no longer symmetric but remained under 
one foot.This is demonstrated at Figure 3-14¢.A crossrange 
dispersion contour was constructed at Figure 3-15 but this 
did not emphasize the resonance phenomenon sufficiently due 
to linearized aerodynamics and small dispersion 
Magnitudes.The parameters that varied most significantly 
near resonance were angle of attack and burn time.Indeed,it 
is clear that the projectile is unstable in this region with 
€a maximum angle of attack of 1.21 radians(7/70 degrees) .Due to 
the increased angle of attack, drag increased and burn time 
decreased in the order of seven seconds (5%) compared to the 
baseline trajectory. This clearly affected the range at 
burnout.Due to large angle of attack the linearized 
aerodynamics are not valid.Projectile performance at wind 
resonance is worthy of additional analysis.It was decided to 
demonstrate resonance by plotting the maximum angle of 
attack Calpha) versus wavelength with different thrust/ drag 
scenarios.The resonance results are at Figure 3-16.The 
projectile natural frequency is estimated at 0.22 Hz. 


Resonance was confirmed in the 0.22 Hz regime for thrust 


errors under 10% including varying Vwo.This was necessary as 
the system is non-linear and Vwo could have affected 
resonance.It was noted that the curve shifted and increased 
or decreased in magnitude slightly in each case.This is due 
to the non-linearity of the system equations and the 
variations in Vo during flight.The absolute case of no 
thrust showed resonance at 10000 feet with a maximum angle 


of attack of 1.25 radians(72 degrees). 
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Figure 3-11 Dispersion Versus Wind Phase Angle(Large Lambda) 
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Figure 3-12 Dispersion Versus Phase Angle 


(Medium Lambda) 
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Figure 3-13 Dispersion Versus Phase Angle(Small Lambda) 
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Figure 3-14 Dispersion Versus Phase Angle 
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Figure 3-15 Phase Angle Versus Wavelength 


3025 
*10° 


WL 
OQ V¥O = 60 | 
- : 
a 
e 
A 
eT 
(Gh 
HOO 
er 
ie i 
® fe ny 
= . ‘ 
ah 
f AW 
} Gy 
TS d : | 
Y he I 
a : : 
ai 
<— Ww fut Y 
mea se A\ 
0° 
<f { 
eal | 
f \ 
> en 
\ 
<— } 
cS ! . 
LO { . 
<j S / a) : 
= ! 
CL, \ ‘ | 
Lp : 
\- 
nd vA 
<q / ae 
/ we 
} \ ‘\ | 
f De | 
‘ \ 
° XN 
> = ry ) 
= <a Ly | 
ee 
SS ee 
oe 
ere = : 
| 
| 
| 





= 
8000 8500 9000 9500 10000 10500 11000 11500 12000 
WAVELENGTH (FEET) 


Figure 3-16 Detrimental Resonance Between Projectile 
Short Period and Spatial Frequency of 


Wind Profile 


It is clear that if the wind profile experienced by 
the projectile has harmonic content near the natural 
frequency then the round will be unstable and/or have large 
dispersion at burnout. This is best illustrated by viewing 
the wind as Vwtt) with reference to the projectile. By a 
Fourier transform of the wind profile resonance with 
projectile natural frequency (fn) will be evident in the 
frequency domain.The wind oscillations need only be velocity 
variations with the appropriate harmonic content.Figure 3- 
17, which illustrates this concept, shows fn coinciding with 
a maximum of Vwtf). One would anticipate stability problems 


for this case. 
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Figure 3-17 Fourier Transform Of Wind Profile Concept 
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Calculations for the long period resonance or 
phugoid show that the phugoid is irrelevant as flight time 
ig too short: 

T=4.44*Vo0/g=301.8 seconds (Reference 3] 

The results of the sinusoidal winds were also 
compared with the aerodynamic jump phenomenom outlined in 
BRL report 1077 CReference 4 J]. At first glance the 
Sinusoidal wind results appear consistent with aerojump,bdut 
detailed calculations show otherwise. The rate of change of 
initial angle of attack is critical to aerojump and is 
easily incorporated into the aerojump equation.The results 
for the sine wind are a function of wavelength and of the 
same form as the simplified model results but the cosine 
wind and constant wind both yielded a small aerojump 
proportional to the initial wind.The magnitudes of the 
aerojump dispersions were more consistent with the net lift 
transient effects. Further, continuous perturbations in the 
form of wind after time zero are not in the aerojump 
equations and for all practical purposes the phenomenom is 
not observed in the model.These results should have been 
apparent before calculations were performed. Aerojump 
equations are not appropriate nee non-zero wind 
conditions.Clearly, when wind is included the aerojump 


equations cannot predict dispersion. 
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B. FIVE DEGRERS) OF (FREEDOM 

As previously stated dispersion was measured from the 
ideal pseudo-vacuum trajectory at a fixed time £07; 
initiation of terminal guidance (Tc) 101 seconds.Thrust 
error sensitivity in the zero wind condition is presented in 
Figures 3-18 and 3-19.In these two figures it is clear that 
thrust errors can have a significant impact on the location 
of ~ the progecti le -atet hice sera gute indicates that as 
thrust increases the dispersion moves downrange as 
expected.Figure 3-19 shows that as thrust increases the 
height of the projectile at Tc decreases. Due to added 
thrust the projectile is located at a point further along 
the parabolic trajectory; hence the height is less. 

As well Figure 3-20 shows these effects on available 
burn time of the ramjet.For thrust equals drag burnout at 
101 seconds exactly equals the Tc.With less thrust, terminal 
guidance is initiated before all the propellant is’ burned. 
Conversely, when thruat exceeds drag, burnout occuréd jon ten 
to initiation of terminal guidance.Assuming this fixed time 
for initiation of terminal guidance, the projectile may 
experience up to seven aeconds of unpowered flight ¢ 7% 
error) before the guidance ig deployed. The period of 
unpowered flight with associated high drag will cause major 


dispersions as shown in subsequent figures. 
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Figure 3-20 Thrust Error Effects (Burn Time) 


Crosswind sensitivity results were very much the same 
as the three degree simulation.In fact the results appear 
even more linear with the same rotation about the zero zero 
point as in Figure 3-21.The gross features of Figure 3-3 and 
3-21 are identical but Figure 3-3 shows dispersion at 
burnout whereas Figure 3-21 shows dispersion at Tc. 

The crossrange dispersion results of a no thrust 
projectile with constant crosswind profiles and with thrust 
equal to drag (transients effects) are at Tables 3.3 and 3.4 


respectively. 


Table 3.3 No Thrust Wind Effects 
Five Degrees 


Thrust =0 time=101 seconds 


Vwxy (feet/sec) dispersion (feet) 
crossrange 
0 0 
30 S20 
60 185.11 
So) 272,90 
120 370797 
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Figure 3-21 Crosswind Effects Five Degrees 


Table 3.4 Dispersion Due to Net Lift 
from Transients 


Thrust=drag time=101 secs 


Vwxy (feet/sec) dispersion(feet) 
O 0 
30 -9202 
60 0413 
a6) .0643 
120 <Osd9 


Burn time was not significantly affected by crosswinds. 
Headwind and tailwind effects on range,height' and 
burntime in different thrust scenarios are shown in Figures 
3-22 to 3-26. It is clear that dispersion from the pseudo- 
vacuum trajectory iS Mminimized when thrust nearly equals 
drag.It is also apparent that the overthrusted trajectories 
involve greater dispersion aS burnout occurS prior to 
terminal guidance and unpowered flight occurs for some time 
dependent on the scenario.This same error is observed in the 
case of headwinds and is most clear in the Figures 3-22B to 
3-25B about the thrust equal to drag axis. (Note that Figure 
3-22 to 3-25 inclusive have A,B and C sub-figures onder are 
Simply enlargements of key portions of the respective main 


figure). 
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Figure 3-22 Headwind Sensitivity (Range) 
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Figure 3-224 Headwind Sensitivity (Range) 
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Figure 3-22B Headwind Sensitivity (Range) 
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Figure 3-22C Headwind Sensitivity (Range) 
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Figure 3-23 Tailwind Sensitivity (Range) 
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Figure 3-23A Tailwind Sensitivity (Range) 
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Figure 3-23B Tailwind Sensitivity (Range) 
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Figure 3-23C Tailwind Sensitivity (Range) 
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APPENDIX D PSEUDO-VACUUM TRAJECTORY PROGRAM 


DIM X(SOO),Y( 3am), RO(Saa),D( saa), MDOT (Sama), mM( 3am), V (saa) 
DIM CDO( 3am), MA(Saa) 
DIM T(3aa) 
MMIN=198 
ISP=95a 
DELT=.5 
M(1) S225 
THETA=. 7852981634007006# 
VO=2452 
LPRINT "TIME", "X","Y", "MASS", "SPEED" 
FOR I=1 TO 30a 
X (1) =VO*COS (THETA) *T (1) 
Y (1) = (VO#SIN (THETA) *#T (1) )— (16. 287*"T (1) *2) 
V (I) =( (VO*#COS (THETA) ) “2+ (VO*#SIN (THETA) — 
Ze. 174*T(1))*2)*.5 
TEMP =518. 688-. aas5662*Y (1) 
RO(I)=(2116. 22* (1-—(. AASS5E62*Y (1) /518. 688) )~S. f56) / 
(1716. 5*TEMP) 
IF Y(I))=36289! THEN RO(I)=. aa1271 
IF Y(1I))=36289! THEN TEMP=389. 99 
MA(I)=V(1)/SQR(1.4%1716. 3* TEMP) 
CDO( 1) =. SM@15—. 1@25*MA(T) 
D(1)=. 276%. S*RO(1) *. 3491*V (1) 2 
D(1I)=D(1) /.276*CDO(I) 
MDOT (I+1)=D(1)/ISP 
M(I+1)=M(1)-MDOT (I) *DELT 
T(I+1)=T(1)+DELT 
arNiied (lex (1.).¥(1),M(1), Vil) 
IF M(I+1) <=MMIN GOTO 222 
NEXT I 
LPRINT "RAMJET BURNOUT BETWEEN T="sT(I) s"AND";T(I+1) 


WwW 
fi 


burning.Burnout will occur on the vacuum trajectory but not 


at the same point, as shown in Figure C-2. 


Vwt 
Vwe 








Figure C-2 Burnout Variations 


el 


In Figure C-1 wind (Vw) has been added for the curve #4 
in the x direction.The new equations are 


mx4 = (-Dy + Dw) + Tx 


mz4 -D>s + Ts -mg 


where Dx + Dw= 1/2 9 Ue ea & 
and Uxtwind relative to the projectile) =Vx+Vw. With no 
thrust curve 4 results. 


Finally, if thrust is added to equal drag the 


following equations are valid: 


mXxS = 0 or Ty = Dx + Dw 
NZS = -mg or Ts = Dg 


The accelerometer(s) read the acceleration caused by any 
differences in the thrust and drag forces (zero in this 
case). 

It is clear that curves 1 and 5 are equal (and 3 if 
thrust=drag) but the burntime of a ramjet will not be equal 


in each case due to different drag profiles. 


t = (mass fuel) ¢Isp)/thrust 


Therefore headwinds will require more thrust and 


burntime will be shorter.Conversely tailwinds allow longer 


31 


case Ux=Vx and Uz=Vz.This trajectory is represented by curve 
2.-An accelometer in the x direction measures Dx/m and Dz2/m 
in the z direceven- 

If thrust is added, the trajectory is as shown at curve 


3(thrust greater than drag).The key equations are then: 


mxX3 = -Dy + Ty 
MZg = -Ds + Ts -MgG 
Accelerometers would measure (Tx-Dx)/m and (Ty-Dy)/m 


respectively.If thrust equals drag the trajectory would be 


curve 1. 





Figure C-1 Pseudo-Vacuum Trajectory Concept Curves 
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APPENDIX C-PSEUDO-VACUUM TRAJECTORY CONCEPT 


The basic concepts behind the pseudo-vacuum are not 
always readily apparent.This appendix analyzes the basic 
equations of motion in the following sequence of cases to 
jlluatrates the concept: 

vacuum trajectory 
atmosphere trajectory 
thrusted trajectory 
wind influenced trajectory 
thrust equal drag trajectory 
In Figure C-1 (curve 1) the applicable equations of motion 


are: 
0 


a 
ms 

= 

i 


J 
N 
~ 
WL 


An accelometer on the projectile in the z direction would 
measure no acceleration as it is in free fall already. This 
is similar to a satellite orbit or someone on an elevator in 
free fall. 


When atmosphere is added the key equations are: 


miv> 


=e = 47 ? Uyx* A Ca 


MZ? 


-Dz - mg 


where U is relative wind and D is atmospheric drag.In this 
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When t=O, the equation reduces to: 

-2§w +1(28w) = O which is expected if the initial 
condition is zero.Cie. dispersion at zero equals zero at 
time zero). 

As t approaches infinity the equation value approaches 
infinite dispersion, but the oscillations of alpha have 
damped out. Note that the first derivative, proportional to 


velocity, tends to a constant value: 


(w/(§2 + we ))/m as t approaches infinity which is 


also to be expected. 


APPENDIX B TRANSIENT ANGLE OF ATTACK CALCULATIONS 


The following equation of motion applies to the dynamics 
of angle of attack (Calpha) oscillations due to the torque 


equation described in the three-degrees-of-freedom model. 


my = Lw which is proportional to alphact) 
where angle of attack is of the form 
-st 
e< (t)=e sin(wt) where 8is the damping constant 


and w is the angular velocity. Therefore: 


7 
) ax<c¢ct) dt dT 
foe °) 


Ae | 
5 § e-St sin(wt) dt aT 
0 ”*0O 


my 





e-ST 






(-Ssin(wT) - w cos(wT)) + 
+ we 


we) dT 





= wt(S 2 +w2) -28w + 


e-St(§2sin(wt) +2 § woos(wt) -w2sin(wt)) 





CIS 20te 2) 2 
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impulse=950 “seconds” 

Mass Property Equations: 
CG= (225x39.4)-(32.0xfuel burned)/(223-fuel burned aan 
Izz=.35-¢.0019xfuel burned) in slug-ft 


Ixx=Iyy=18.8-¢(.0423xfuel burned) in slug-ft 


Aerodynamic Coefficient Equations: 
Cd=Cdo +.5 x Cl x sin &) 
Cdo= —-.1025 x NACH e+e. >0l> Cimaeo< MACHT < 222 
Cm=Cl (Xcg-Xcp)/d 


Ci=Cl x ad 


Key Wind Tunnel Results (per degree): 
MACH=1.8 


Cdo=.317 Cmq=-203 


ox Oo Z ef 6 8 10 
GH 0 aS oo 62 232 Ls 
ACD 45.4 45.6 48.2 49.7 30.2 350.0 
MACH=2.2 
Cdo=.276 Cmq=-190 
x 0 2 4 6 8 10 
Ca 0 -19. .4009Feo aca as 


Xcp 42.9 42.9 44.6 45.1 46.0 44.6 
Cmqzd(Cm)/dt(qd/2Vv) 


Initial Muzzle Spin Rate 30 PI radians /sec 


APPENDIX A-PROJECTILE CHARACTERISTICS 
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Figure A-1 ARMT Projectile Baseline Design 


The following characteristics are obtained from 


References 2 and 6. 

length=72 inches 

diameter=8 inches 

launch weight=225 lbm 

launch CG=39.4 inches from spike apex 


fuel weight=27 lbm (CG at 52.0 inches from apex) 
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except where harmonic content near resonance is 
contained in the wind profile. 

-headwinds create greater dispersion as ramjet 
burnout occurs before terminal guidance is 
initiated.Crosswinds and tailwinds generally will 


have minimal effects. 


Recommendations for further study related to this 
thesis are: 
-Dynamics at launch including the initial spin 
damping, Magnus effects and aerodynamic jump. 
-Dynamics at terminal guidance including air inlet 
spike ejection and canard deployment. 
-Complete adaptation of BRL HTRAJ on the NPGS ' IBM 


system to use non-linear aerodynamics. 


G3 
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IV. CONCLUSIONS 

Generally this thesis confirms the viability of the 
pseudo-vacuum trajectory from a thrust equals drag 
projectile.Indeed winds up to 120 feet per second and thrust 
errors up to five percent created dispersions of less than 
600 feet over a range of 175090 feet for the ARMT 
projectile. 

Specifically the following conclusions are made: 

-the drag due to wind (Dw) is the aerodynamic 
force which creates the most dispersion. 

-the net lift due to transients (Lw) is a minor 
aerodynamic force in creating dispersion. 

-the dispersion is proportional to the damping 
coefficient (Cmq),although for practical values of Cmq 
dispersion can be treated as insensitive to Cmq. 

-the projectile resonance or short period is at 
approximately OO. 22 hz and is unstable in the 
resonance regime due to variable wind profiles.The phugoid 
is irrelevant. 

-thrust errors are the most significant errors 
contributing to dispersion. 

-winds have a minor effect on dispersion during 


powered 23 IPE p ove, when thrust approximates drag, 


Co 
fit 





Figure 3-28 Three Dimensional Trajectory With Altitude 
Dependent Wind Profile 
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Wind effects during near thrust equals drag portions of 
the trajectory were very small.The erlie2zeal fect Gre ais 
affecting the dispersion will be the setting of the terminal 
guidance clock in the projectile.This presents’ several 
options for ‘setting initiation of terminal guidance: 

-Fixed time irrespective of trajectory Cie. 
elevation of gun).Burn time would be shorter at less than 45 
degree launch and the basket would be more sensitive to wind 
and thrust errors.As well the trajectory would change and so 
would the predicted basket. 

-Fixed time for a given trajectory.This requires 
analysis of the sensitivity to elevation of the gun and a 
manual Tc setting method for gunnersg. 

-Ramjyet burnout.This would ensure that Tc occurred 
at an “optimum™ point almost exactly on the pseudo-vacuum 
trajectory but the location along it at Tc would be 
Per icult to predict. 

-At given pressure for altitude (Reference 5 ].This 
would establish accurate baskets for given gun elevations 
assuming the altitude chosen was before ramjet burnout. 

several trajectories were Simulated with a varying wind 
profile aS a function of altitude in the Standard 
Atmosphere.These are shown at Figure 3-28 as compared to the 


pseudo-vacuum trajectory. 
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Velocities at time 101 seconds were measured and are 
shown at Figure 3-27.The total variation of extreme cases is 


from 2285 feet/sec to 2315 feet/sec, or only 30 feet/sec. 
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Figure 3-27 Velocity Sensitivity 
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Deviations from the ideal trajectory were correlated 
and tabulated in Table 3.5.ResultsS were conservatively 
linearized and are presented in that manner. Wind errors are 
read in addition to the thrust errors at the left of the 
table.For example: 


If thrust error=-5%* 
crosswind=120 ft/sec 


predicted dispersion is X=-540 ft 


Vetere 
Z=130 ft :from ideal basket. 


Table 3.5 Dispersion Analysis Results 


Thrust Errors Winds 
continuous 120 ft/sec 


Cross Head Tail 
Errors »,4 Zz Y xX iz X Z 


QO % oe 1 
1 % 10 1 
2 % 20 ez 
> % =] 6) 10 
10 % 100 30 





ie 


TAIL/HEADWIND VELOCITY(FT/SEC) 


—60 


120 


90 


60 


30 


0 





-30 


LEGEND 
THRUST=DRAG 


90 


meme eee eee eee 


TIME=101 SECS 


-—120 


85 90 95 100 105 110 
BURN TIME (SECONDS) 


Figure 3-26 Tail/Headwind Effects (Burntime) 
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Figure 3-25C Tailwind Sensitivity (Height) 
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Figure 3-25B Tailwind Sensitivity (Height) 
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Figure 3-25A Tailwind Sensitivity 
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Figure 3-25 Tailwind Sensitivity (Heron. 
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Figure 3-24C Headwind Sensitivity (Height) 
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Figure 3~24B Headwind Sensitivity 
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Figure 3-24A Headwind Sensitivity (Height) 
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Figure 3-24 Headwind Sensitivity (Height) 
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